Ancient mitogenomes from Pre-Pottery Neolithic Central Anatolia and the effects of a Late Neolithic bottleneck in sheep (Ovis aries).

Occupied between ~10,300 and 9300 years ago, the Pre-Pottery Neolithic site of Asikli Höyük in Central Anatolia went through early phases of sheep domestication. Analysis of 629 mitochondrial genomes from this and numerous sites in Anatolia, southwest Asia, Europe, and Africa produced a phylogenetic tree with excessive coalescences (nodes) around the Neolithic, a potential signature of a domestication bottleneck. This is consistent with archeological evidence of sheep management at Asikli Höyük which transitioned from residential stabling to open pasturing over a millennium of site occupation. However, unexpectedly, we detected high genetic diversity throughout Asikli Höyük's occupation rather than a bottleneck. Instead, we detected a tenfold demographic bottleneck later in the Neolithic, which caused the fixation of mitochondrial haplogroup B in southwestern Anatolia. The mitochondrial genetic makeup that emerged was carried from the core region of early Neolithic sheep management into Europe and dominates the matrilineal diversity of both its ancient and the billion-strong modern sheep populations.

Assessing the performance of Moghani crossbred lambs derived from different mating systems with Texel and Booroola sheep.

In our ongoing project, which focuses on the introgression of Booroola/FecB gene and the myostatin (MSTN) gene into purebred Moghani sheep, we assessed the performance of second-generation Moghani crossbreds such as second crossbreds (F2) and initial backcross generation (BC1). These crossbreds were generated through different mating systems, including in-breeding, outcrossing, first paternal backcrossing (PBC1), and first maternal backcrossing (MBC1). Notably, F2 strains exhibited lean tail, woolly fleece and a higher percentage of white coat color compared to BC1. The impact of mating systems and birth types on pre-weaning survival rates was found to be statistically significant (P < 0.0001), with singleton offspring resulting from paternal backcross showing a particularly substantial effect. The F2 crossbred lambs carrying the Booroola gene did not show a statistically significant difference in survivability compared to those carrying the MSTN gene, implying the Booroola prolificacy gene had no significant impact on survival outcomes. However, the occurrence of multiple births had a significant negative impact on lamb survival (P < 0.0001). The PBC1 sheep strains, specifically Texel Tamlet ram strains carrying the MSTN mutation, exhibited superior growth rates compared to others (P < 0.05). Interestingly, the MSTN mutation in the homozygous variant genotype significantly impacts growth rate before weaning compared to other genotypes and pure Moghani sheep (P < 0.05). In conclusion, this study objectively underscores the pivotal role of genetic factors, specifically through strategic mating systems like paternal backcrossing, in enhancing desired traits and growth rates in Moghani sheep, thereby contributing valuable insights to the field of sheep breeding programs.

Screening of hair follicle telogen-associated circRNAs in sheep and construction of their ceRNA network.

Sheep breeds with hair-shedding traits have many advantages over non-shedding sheep breeds, not only because of reduced shearing labor and feeding management costs but also because it reduces in vitro parasites and improves adaptability to summer heat stress. The wool of Dorper sheep naturally sheds in spring due to the periodic growth of hair follicles. CircRNAs primarily regulate the morphogenesis of hair follicles through the ceRNA mechanism. In this study, five 2-year-old Dorper ewes with extreme hair-shedding phenotype (S) and three Dorper ewes with non-shedding (N) phenotype were selected for subsequent analyses. For RNA extraction, skin tissues were collected on 27th September 2019 (S1, N1), 3rd January 2020 (S2, N2), and 17th March 2020 (S3, N3), which were then subjected to RNA-seq. RNA-seq technology revealed 20,185 novel circRNAs in the hair follicles of Dorper sheep. Among them, 1450 circRNAs were differentially expressed (DE). Clustering heatmap and expression pattern analyses were performed on DE circRNAs, which indicated 78 circRNAs with T pattern (Telogen, highly expressed in telogen), and the source genes for candidate circRNAs were further screened by functional enrichment analysis, which identified 13 crucial genes enriched in pathways associated with hair follicle development. Additionally, a ceRNA regulatory network comprising 4 circRNAs, 11 miRNAs, and 13 target genes was constructed. Overall, this study screened circRNAs that may be associated with the telogen phase of hair follicles in sheep, providing a relevant theoretical basis for wool shedding in sheep and for breeding Dorper sheep with automatic wool shedding.

Inheritance of breeding traits in Karakul sheep under different selection approaches.

In the previous century, each intrabreed type of the Karakul sheep breed was characterized by significant numbers, representing a super population with rich genetic diversity. However, over time, the genetic diversity within the breed's gene pool has undergone significant depletion. At present, the Karakul breed is predominantly composed of only two small populations, distinguished by their fur colors: black and gray. Consequently, under such circumstances, genetic advancements in breeding endeavors are likely to be relatively limited, especially given the potential risk of these populations disappearing altogether in the future. Hence, the preservation and judicious utilization of the available genetic resources within the black and gray Karakul sheep populations hold paramount importance in breeding efforts. The primary objective of our research was to investigate the heritability of breeding traits among gray lambs through various selection options. The study was conducted at the "Kumkent" base farm in the Sozak district of the Kyzylorda region. Our findings revealed that the inheritance of gray and black fur colors across the different selection options occurred in a consistent ratio. In the first selection variant (a gray ram with even silver marking ♂ x a black jacket fur type ewe with intense pigmentation ♀), the proportion of gray offspring was 50.6%. Similarly, in the second selection variant (a black jacket fur type ram with intense pigmentation ♂ x a gray ewe with even silver marking ♀), the proportion of gray offspring was 49.6%. The percentage of black lambs obtained in both selection options was nearly equivalent, with 49.4% and 50.4% in the first and second variants, respectively.

Identification and screening of circular RNAs during adipogenic differentiation of ovine preadipocyte by RNA-seq.

Circular RNAs (circRNAs) are a class of non-coding RNAs that play important roles in preadipocyte differentiation and adipogenesis. However, little is known about genome-wide identification, expression profile, and function of circRNAs in sheep. To investigate the role of circRNAs during ovine adipogenic differentiation, the subcutaneous adipose tissue of Tibetan rams was collected in June 2022. Subsequently, the preadipocytes were immediately isolated from collected adipose tissue and then induced to begin differentiation. The adipocytes samples cultured on days 0, 2, and 8 of preadipocytes differentiation were used to perform RNA sequencing (RNA-seq) analysis to construct the expression profiles of circRNAs. Subsequently, the function of differentially expressed circRNAs was investigated by performing the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of their parent genes. Finally, a circRNAs-miRNAs-mRNAs network involved in adipogenic differentiation was been analyzed. As a result, a total of 6,449 candidate circRNAs were identified in ovine preadipocytes. Of these circRNAs identified, 63 candidate circRNAs were differentially expressed among the three differentiation stages and their parent genes were mainly enriched in acetyl-CoA metabolic process, positive regulation of lipid biosynthetic process, positive regulation of steroid biosynthetic process, and focal adhesion pathway (P < 0.05). Based on a circRNAs-miRNAs-mRNAs regulatory network constructed, circ_004977, circ_006132 and circ_003788 were found to function as competing endogenous RNAs (ceRNAs) to regulate ovine preadipocyte differentiation and lipid metabolism. The results provide an improved understanding of functions and molecular mechanisms of circRNAs underlying ovine adipogenesis in sheep.

The moderate fat deposition contributes to improve mutton quality, which is associated with the differentiation of preadipocytes. To investigate roles of circular RNAs (circRNAs) in preadipocyte differentiation, we identified circRNAs on days 0, 2, and 8 of preadipocytes differentiation and compared the expression profile of circRNAs at different adipogenic differentiation stages. A total of 6,449 candidate circRNAs were identified, among which 63 candidate circRNAs were differentially expressed among the three differentiation stages. The parent genes of differentially expressed circRNAs were enriched in several biological process and pathways related to lipid metabolism and synthesis. In addition, several circRNAs may regulate ovine preadipocyte differentiation by interacting with microRNAs (miRNAs). The results reveal the potential roles of circRNAs in adipogenic differentiation of sheep.

Increasing the meat productivity of young sheep based on the use of the gene pool of the Dorper and Hissar breeds.

In the pursuit of enhanced mutton production, improving the genetic reservoir of sheep with early maturation and high meat productivity is imperative. This study aims to assess the efficacy of integrating Dorper and Hissar rams into the breeding program of Kazakh fat-tailed coarse-haired ewes for generating young mutton. The research involved forming three groups, each comprising 40 ewes of the Kazakh fat-tailed coarse-haired breed, based on analog pairs. Ewes in Group I were inseminated with Dorper ram semen, those in Group II were inseminated with Hissar ram semen, and Group III served as a control with purebred Kazakh fat-tailed coarse-haired sheep breeding. Results revealed that crossbred rams in Group II achieved a significantly higher live weight of 45.2 kg at 120 days of age, surpassing the other groups by 9.7 kg and 10.6 kg. Crossbred gimmers in Group II reached a live weight of 42.0 kg by 4 months, outpacing the other groups by 12.2 kg. The crossbred lambs exhibited an expansive, deep, and sturdy physique, indicative of elevated meat productivity. Physique index analysis displayed that crossbred rams exhibited elongated limbs, bulkiness, and massiveness compared to purebred Kazakh fat-tailed coarse-haired lambs. In the 4.0-4.5-month age range, crossbred rams demonstrated a higher carcass muscle yield than their purebred counterparts, albeit the latter exhibited a 0.18% greater bone yield. Moreover, the meat of groups I and II sheep contained 19.6% and 20.1% protein content, respectively, surpassing the local Kazakh fat-tailed sheep population by 0.7% and 1.2% in absolute terms.

Investigating pedigree- and SNP-associated components of heritability in a wild population of Soay sheep.

Estimates of narrow sense heritability derived from genomic data that contain related individuals may be biased due to the within-family effects such as dominance, epistasis and common environmental factors. However, for many wild populations, removal of related individuals from the data would result in small sample sizes. In 2013, Zaitlen et al. proposed a method to estimate heritability in populations that include close relatives by simultaneously fitting an identity-by-state (IBS) genomic relatedness matrix (GRM) and an identity-by-descent (IBD) GRM. The IBD GRM is identical to the IBS GRM, except relatedness estimates below a specified threshold are set to 0. We applied this method to a sample of 8557 wild Soay sheep from St. Kilda, with genotypic information for 419,281 single nucleotide polymorphisms. We aimed to see how this method would partition heritability into population-level (IBS) and family-associated (IBD) variance for a range of genetic architectures, and so we focused on a mixture of polygenic and monogenic traits. We also implemented a variant of the model in which the IBD GRM was replaced by a GRM constructed from SNPs with low minor allele frequency to examine whether any additive genetic variance is captured by rare alleles. Whilst the inclusion of the IBD GRM did not significantly improve the fit of the model for the monogenic traits, it improved the fit for some of the polygenic traits, suggesting that dominance, epistasis and/or common environment not already captured by the non-genetic random effects fitted in our models may influence these traits.

Genome-wide association studies of parasite resistance, productivity and immunology traits in Scottish Blackface sheep.

Gastrointestinal parasitism represents a global problem for grazing ruminants, which can be addressed sustainably by breeding animals to be more resistant against infection by parasites. The aim of this study was to assess the genetic architecture underlying traits associated with gastrointestinal parasite resistance, immunological profile and production in meat sheep, and identify and characterise candidate genes affecting these traits. Data on gastrointestinal parasite infection (faecal egg counts for Strongyles (FECS) and Nematodirus (FECN) and faecal oocyst counts for Coccidia, along with faecal soiling scores (DAG), characterised by the accumulation of faeces around the perineum) and production (live weight (LWT)) were gathered from a flock Scottish Blackface lambs at three and four months of age. Data on the immune profile were also collected from a subset of these lambs at two and five months of age. Immune traits included the production of Interferon-γ (IFN-γ), Interleukin (IL)-4 and IL-10 following stimulation of whole blood with pokeweed mitogen (PWM) or antigen from the gastric parasite Teladorsagia circumcincta (T-ci), and serum levels of T. circumcincta-specific immunoglobulin A (IgA). Animals were genotyped with genome-wide DNA arrays, and a total of 1 766 animals and 45 827 Single Nucleotide Polymorphisms (SNPs) were retained following quality control and imputation. Genome-wide association studies were performed for 24 traits. The effects of individual markers with significant effects were estimated, and the genotypic effect solutions were used to estimate additive and dominance effects, and the proportion of additive genetic variance attributed to each SNP locus. A total of 15 SNPs were associated at least at a suggestive level with FECS, FECN, DAG, IgA, PWM-induced IFN-γ and IL-4, and T-ci-induced IL-10. This study uncovered 52 genes closely related to immune function in proximity to these SNPs. A number of genes encoding C-type lectins and killer cell lectin-like family members were close to a SNP associated with FECN, while several genes encoding IL-1 cytokine family members were found to be associated with IgA. Potential candidate genes belonging to or in close proximity with the Major Histocompatibility Complex (MHC) were revealed, including Homeostatic Iron Regulator and butyrophilin coding genes associated with IFN-γ(PWM), and IL-17 coding genes associated with IgA. Due to the importance of the MHC in the control of immune responses, these genes may play an important role in resistance to parasitic infections. Our results reveal a largely complex and polygenic genetic profile of the studied traits in this Scottish Blackface sheep population.

Wool production and quality traits of pure- and crossbred Merino-type sheep.

Very little information is available on the quality of wool produced by terminal crosses out of wool producing dam lines. This study was therefore undertaken to elucidate the wool quality traits of four crossbred lines relative to Merinos and Dohne Merinos, which served as dam lines. Each dam line was mated to rams of their own breed as well as to Dormer or Ile de France rams to create four crossbred lines, namely, Dohne x Dormer, Dohne x Ile de France, Merino x Dormer, and Merino x Ile de France, in addition to the purebred Dohne Merino and Merino lines. Four rams and four ewes from each of these lines were reared up to one year of age under optimal growth conditions and shorn at the end of that time. Fleeces were weighed and samples collected for analysis. Neither sex nor genotype influenced clean fleece weight or clean yield percentage, but fibre diameter differed between genotypes. Purebred Merino had the finest wool (18.26 µm) and the Merino x Dormer cross the coarsest (26.01 µm). However, all lines still produced wool that could be used in manufacturing clothing, while fleeces showed good uniformity. The two purebred lines had the highest comfort factor (> 98%). The crossbred lines produced a similar quantity of wool as the purebreds, although of poorer quality. All genotypes except Merino x Dormer appear to produce wool that could be processed into garments, although the crossbred wool would only be suitable for outerwear.

Phenotype transition from wild mouflon to domestic sheep.

The domestication of animals started around 12,000 years ago in the Near East region. This "endless process" is characterized by the gradual accumulation of changes that progressively marked the genetic, phenotypic and physiological differences between wild and domesticated species. The main distinctive phenotypic characteristics are not all directly attributable to the human-mediated selection of more productive traits. In the last decades, two main hypotheses have been proposed to clarify the emergence of such a set of phenotypic traits across a variety of domestic species. The first hypothesis relates the phenotype of the domesticated species to an altered thyroid hormone-based signaling, whereas the second one relates it to changes in the neural crest cells induced by selection of animals for tameness. These two hypotheses are not necessarily mutually exclusive since they may have contributed differently to the process over time and space. The adaptation model induced by domestication can be adopted to clarify some aspects (that are still controversial and debated) of the long-term evolutionary process leading from the wild Neolithic mouflon to the current domestic sheep. Indeed, sheep are among the earliest animals to have been domesticated by humans, around 12,000 years ago, and since then, they have represented a crucial resource in human history. The aim of this review is to shed light on the molecular mechanisms and the specific genomic variants that underlie the phenotypic variability between sheep and mouflon. In this regard, we carried out a critical review of the most recent studies on the molecular mechanisms that are most accredited to be responsible for coat color and phenotype, tail size and presence of horns. We also highlight that, in such a complicate context, sheep/mouflon hybrids represent a powerful and innovative model for studying the mechanism by which the phenotypic traits related to the phenotypic responses to domestication are inherited. Knowledge of these mechanisms could have a significant impact on the selection of more productive breeds. In fact, as in a journey back in time of animal domestication, the genetic traits of today's domestic species are being progressively and deliberately shaped according to human needs, in a direction opposite to that followed during domestication.

Ruminal microbiota-host crosstalks promote ruminal epithelial development in neonatal lambs with alfalfa hay introduction.

Ruminal microbiota is gradually established after birth, while microbiota maturation could be highly diverse because of varied solid dietary accessibility. However, how the ruminal microbiota accreted from postnatal hay diets alters rumen epithelial development, and how this affects animal health remains largely unknown. Here, neonatal lambs were introduced to starchy corn-soybean starter or corn-soybean starter + alfalfa hay (AH) to investigate the influences of early life ruminal microbiome on rumen epithelial development using integrated 16s rRNA sequencing-metagenome-transcriptome approaches. The results showed that AH introduction elevated average daily weight gain, rumen weight and volume, rumen epithelial papillae length, and rumen muscle layer thickness. Meanwhile, the relative abundance of fibrolytic bacteria (Christensenellaceae R-7 group, Prevotellaceae UCG-001, and Succinivibrio), acetate producer (Acetitomaculum and Mitsuokella), and propionate producer Succiniclasticum was increased in the rumen content by AH supplementation (P < 0.05). Moreover, AH introduction decreased the relative abundance of total CAZymes, CBM, and GH and increased the abundance of KO genes related to volatile fatty acid (VFA) generation in the rumen content. AH lambs had a higher relative abundance of Succiniclasticum, Megasphaera, Succinivibrio, and Suttonella (P < 0.05), while a lower relative abundance of Cloacibacillus, Desulfovibrio, Dialister, Intestinimonas, Parabacteroides, and Pseudoscardovia (P < 0.05) in the rumen epithelial samples. Furthermore, these alterations in ruminal microbial structure and function resulted in ruminal epithelial cell proliferation and development pathways activation. In summary, AH introduction benefited ruminal fiber degradation and VFA generation bacteria colonization and promoted ruminal epithelial development. These findings provide new insights into ruminal microbial-host interactions in the early life.IMPORTANCEWhile it is established that a fiber-rich diet promotes rumen development in lambs, further research is needed to investigate the precise response of rumen microbiota and epithelium to high-quality alfalfa hay. Here, we observed that the inclusion of alfalfa hay led to a discernible alteration in the developmental trajectory of the rumen. Notably, there was a favorable shift in the rumen's volume, morphology, and the development of rumen papillae. Furthermore, ruminal microbial structure and function resulted in ruminal epithelial cell proliferation and development pathways activation, collectively provide compelling evidence supporting the capacity of alfalfa hay to enhance rumen development and health through ruminal micrbiota-host crosstalks. Our findings elucidate the functional response of the rumen to alfalfa hay introduction, providing new insights into strategies for promoting healthy development of the rumen in young ruminants.

Genetic profile of adaptive immune traits and relationships with parasite resistance and productivity in Scottish Blackface sheep.

Gastrointestinal (GI) parasites cause significant production losses in grazing ruminants which can be mitigated by breeding animals resistant to disease. Lymphocyte cytokine production and parasite-specific Immunoglobulin A (IgA) are adaptive immune traits associated with immunity to GI parasites. To explore the utility of these traits for selective breeding purposes, this study estimated the genetic parameters of the immune traits in sheep and assessed their relationship with disease and productivity traits. Whole blood stimulation assays were performed on 1 040 Scottish Blackface lambs at two months of age in 2016-2017. Blood was stimulated with either pokeweed mitogen (PWM), a non-specific activator of lymphocytes, and Teladorsagia circumcincta (T-ci) larval antigen to activate parasite-specific T lymphocytes. The type of adaptive immune response was determined by quantifying production of cytokines interferon-gamma (IFN-γ), interleukin (IL)-4, and IL-10, which relate to T-helper type (Th) 1, Th2 and regulatory T cell responses, respectively. Serum T-ci specific IgA was also quantified. Heritabilities were estimated for each immune trait by univariate analyses. Genetic and phenotypic correlations were estimated between different immune traits, and between immune traits vs. disease and productivity traits that were recorded at three months of age. Disease phenotypes were expressed as faecal egg counts (FEC) of nematode parasites (Strongyles and Nematodirus), faecal oocyst counts (FOC) of coccidian parasites, and faecal soiling score; production was measured as lamb live weight. Significant genetic variation was observed in all immune response traits. Heritabilities of cytokine production varied from low (0.14 ± 0.06) to very high (0.77 ± 0.09) and were always significantly greater than zero (P < 0.05). IgA heritability was found to be moderate (0.41 ± 0.09). Negative associations previously identified between IFN-γ production and FOC, and IL-4 production and strongyle FEC, were not evident in this study, potentially due to the time-lag between immune and parasitology measures. Instead, a positive genetic correlation was found between FOC and PWM-induced IFN-γ production, while a negative genetic correlation was found between FOC and T-ci induced IL-10. Live weight was negatively genetically correlated with IFN-γ responses. Overall, IFN-γ and IL-4 responses were positively correlated, providing little evidence of cross-regulation of Th1 and Th2 immunity within individual sheep. Furthermore, T-ci specific IgA was highly positively correlated with PWM-induced IL-10, indicating a possible role for this cytokine in IgA production. Our results suggest that while genetic selection for adaptive immune response traits is possible and may be beneficial for parasite control, selection of high IFN-γ responsiveness may negatively affect productivity.

Exploring genetic variation in potential indicators of resilience in sheep using fibre diameter measured along the wool staple.

Production animals are increasingly exposed to a wide variety of disturbances that can compromise their productivity, health and well-being. As a result, there is a growing need to be able to select animals that are more resilient to environmental disturbances. Fibre diameter variation measured along a wool staple is expected to contain information about how resilient sheep are to the disturbances of their internal and external environment. This study aimed to develop potential resilience indicators from fibre diameter variation, estimate their genetic parameters and assess whether these traits are genetically correlated across three age stages. The study used 6 140 Merino sheep from the Sheep Cooperative Research Centre Information Nucleus Flocks recorded at yearling, 2 years old, and adult ages. Eight potential traits were defined based on theory, literature and exploratory analysis, which were suggested to capture the animal's ability to resist, respond and recover from potential disturbances. Genetic evaluation of the traits was conducted using pedigree-based animal models. The traits were shown to be low to moderately heritable (0.01-0.33) when examined at each of the three age stages. The potential indicators were generally well correlated with one another within age stages. Further, the genetic correlation between the same trait measured at different age stages was moderate to high between yearling and 2 years old (0.35-0.94) and between 2 years old and adults (0.18-0.70), while slightly lower between yearling and adult estimates (0.09-0.62). These results suggest that selection for resilience indicators from fibre diameter is possible; however, further studies are warranted to refine the trait definitions and validate these indicators against other measures of health, fitness and productive performance.

Distinguishing host responses, extensive viral dissemination and long-term viral RNA persistence in domestic sheep experimentally infected with Crimean-Congo haemorrhagic fever virus Kosovo Hoti.

Crimean-Congo haemorrhagic fever orthonairovirus (CCHFV) is a tick-borne, risk group 4 pathogen that often causes a severe haemorrhagic disease in humans (CCHF) with high case fatality rates. The virus is believed to be maintained in a tick-vertebrate-tick ecological cycle involving numerous wild and domestic animal species; however the biology of CCHFV infection in these animals remains poorly understood. Here, we experimentally infect domestic sheep with CCHFV Kosovo Hoti, a clinical isolate representing high pathogenicity to humans and increasingly utilized in current research. In the absence of prominent clinical signs, the infection leads to an acute viremia and coinciding viral shedding, fever and markers for potential impairment in liver and kidney functions. A number of host responses distinguish the subclinical infection in sheep versus fatal infection in humans. These include an early reduction of neutrophil recruitment and its chemoattractant, IL-8, in the blood stream of infected sheep, whereas neutrophil infiltration and elevated IL-8 are features of fatal CCHFV infections reported in immunodeficient mice and humans. Several inflammatory cytokines that correlate with poor disease outcomes in humans and have potential to cause vascular dysfunction, a primary hallmark of severe CCHF, are down-regulated or restricted from increasing in sheep. Of particular interest, the detection of CCHFV RNA (including full-length genome) in a variety of sheep tissues long after the acute phase of infection indicates a widespread viral dissemination in the host and suggests a potentially long-term persisting impact of CCHFV infection. These findings reveal previously unrecognized aspects of CCHFV biology in animals.

MFSD2A frameshift variant in Kerry Hill sheep with microcephaly.

Microcephaly is a rare neurodevelopmental disorder characterized by reduced skull circumference and brain volume that occurs sporadically in farm animals. We investigated an early-onset neurodegenerative disorder observed in seven lambs of purebred Kerry Hill sheep. Clinical signs included inability to stand or severe ataxia, convulsions, and early death. Diagnostic imaging and brain necropsy confirmed microcephaly. The pedigree of the lambs suggested monogenic autosomal recessive inheritance. We sequenced the genome of one affected lamb, and comparison with 115 control genomes revealed a single private protein-changing variant. This frameshift variant, MFSD2A: c.285dupA, p.(Asp96fs*9), represents a 1-bp duplication predicted to truncate 80% of the open reading frame. MFSD2A is a transmembrane protein that is essential for maintaining blood-brain barrier homeostasis and plays a key role in regulating brain lipogenesis. Human MFSD2A pathogenic variants are associated with a neurodevelopmental disorder with progressive microcephaly, spasticity, and brain imaging abnormalities (NEDMISBA, OMIM 616486). Here we present evidence for the occurrence of a recessively inherited form of microcephaly in sheep due to a loss-of-function variant in MFSD2A (OMIA 002371-9940). To the best of our knowledge, this is the first report of a spontaneous MFSD2A variant in domestic animals.

Amplification of Ovis aries papillomavirus type 2 DNA from an ovine cutaneous fibropapilloma.

Seven of 60 Perendale sheep within a flock developed single or multiple exophytic masses on their distal hind limbs. A mass was excised from one sheep and histological evaluation revealed epidermal and mesenchymal proliferation, papillomavirus-induced keratinocyte changes and marked keratohyalin clumping. Ovis aries papillomavirus type 2 DNA sequences were amplified using PCR.

Sept des 60 moutons Perendale d'un troupeau ont développé des masses exophytiques uniques ou multiples sur leurs membres postérieurs distaux. Une masse a été excisée sur un mouton et l'évaluation histologique a révélé une prolifération épidermique et mésenchymateuse, des modifications kératinocytaires induites par le papillomavirus et une agglutination marquée de kératohyaline. Les séquences d'ADN du papillomavirus Ovis aries de type 2 ont été amplifiées par PCR.

Sete de 60 ovelhas Perendale de um rebanho desenvolveram massas exofíticas na porção distal dos seus membros posteriores. Uma massa foi removida de uma ovelha e a avaliação histopatológica revelou proliferação mesenquimal e epidérmica, alterações queratinocíticas induzidas por papilomavírus e aglomeração queratohialina. Sequências de papilomavírus Ovis aries tipo 2 foram amplificadas utilizando PCR.

Siete de 60 ovejas Perendale dentro de un rebaño desarrollaron masas exofíticas únicas o múltiples en sus extremidades traseras distales. Se extirpó una masa de una oveja y la evaluación histológica reveló proliferación epidérmica y mesenquimal, cambios de queratinocitos inducidos por el virus del papiloma y marcada acumulación de queratohialina. Mediante PCR se amplificaron secuencias de DNA del virus del papiloma Ovis aries tipo 2.

Combined purebred and crossbred genetic evaluation of Columbia, Suffolk, and crossbred lamb birth and weaning weights: systematic effects and heterogeneous variances.

Despite the benefits of crossbreeding on animal performance, genetic evaluation of sheep in the U.S. does not directly incorporate records from crossbred lambs. Crossbred animals may be raised in different environments as compared to purebreds. Systemic factors such as age of dam and birth and rearing type may, therefore, affect purebred and crossbred performance differently. Furthermore, crossbred performance may benefit from heterozygosity, and genetic and environmental variances may be heterogeneous in different breeds and their crosses. Such issues must be accounted for in a combined (purebred and crossbred) genetic evaluation. The objectives of this study were to i) determine the effect of dam age and birth type on birth weight, and dam age and birth-rearing type on weaning weight, in purebred and crossbred lambs, ii) test for heterogeneous genetic and environmental variances in those weights, and iii) assess the impact of including weights on crossbred progeny on sire estimated breeding values (EBV). Performance records were available on purebred Columbia and Suffolk lambs. Crossbred information was available on lambs sired by Suffolk, Columbia or Texel rams mated to Columbia, Suffolk, or crossbred ewes. A multiple-trait animal model was fitted in which weights from Columbia, Suffolk, or crossbred lambs were considered different traits. At birth, there were 4,160, 2,356, and 5,273 Columbia, Suffolk, and crossbred records, respectively, with means (SD) of 5.14 (1.04), 5.32 (1.14), and 5.43 (1.23) kg, respectively. At weaning, on average at 122 (12) d, there were 2,557, 980, and 3,876 Columbia, Suffolk, and crossbred records, respectively, with corresponding means of 39.8 (7.2), 40.3 (7.9), and 39.6 (8.0) kg. Dam age had a large positive effect on birth and weaning weight in pure and crossbred lambs. At birth, however, the predicted effect was larger in crossbred and Suffolk lambs. While an increase in a number of lambs born and reared had a strong and negative influence on birth and weaning weight, the size of the effect did not differ across-breed types. Environmental variances were similar at birth and weaning, but additive variances differed among breed types for both weights. Combining purebred and crossbred information in the evaluation not only improved predictions of genetic merit in purebred sires but also allowed for direct comparisons of sires of different breeds. Breeders thus can benefit from an additional tool for making selection decisions.

Combining multiple breeds in a genetic evaluation allows for their direct comparison. However, differences in management and other systematic effects among breeds may affect the evaluation. Estimates of genetic merit of sires may also be biased by heterosis in crossbred progeny. We examined genetic and environmental factors that affect the efficacy of a multi-breed genetic evaluation. Birth and weaning weights of Columbia, Suffolk, and their cross, were available. Depending on the breed type, the systematic effects of dam age and either birth or birth-rearing type on weights differed. Separately for birth and weaning, weights were defined as a different trait for each breed type. A multi-breed, multi-trait model was fitted that accounted for systematic effects unique to a breed type, and heterosis. Estimated direct and maternal heritabilities were moderate. Genetic correlations between breeds were moderate to high. Estimates of genetic merit of Columbia and Suffolk sires were unaffected by bias due to heterosis and environmental effects when crossbred lambs were included in a purebred or a combined Columbia, Suffolk, and crossbred evaluation. For direct across-breed comparisons, breed type-specific adjustments for systematic effects are necessary when combining weight data on pure and crossbred lambs in a joint genetic evaluation.

Expression profiles of oviductal mRNAs and lncRNAs in the follicular phase and luteal phase of sheep (Ovis aries) with 2 fecundity gene (FecB) genotypes.

FecB (also known as BMPR1B) is a crucial gene in sheep reproduction, which has a mutation (A746G) that was found to increase the ovulation rate and litter size. The FecB mutation is associated with reproductive endocrinology, such mutation can control external estrous characteristics and affect follicle-stimulating hormone during the estrous cycle. Previous researches showed that the FecB mutation can regulate the transcriptomic profiles in the reproductive-related tissues including hypothalamus, pituitary, and ovary during the estrous cycle of small-tailed Han (STH) sheep. However, little research has been reported on the correlation between FecB mutation and the estrous cycle in STH sheep oviduct. To investigate the coding and noncoding transcriptomic profiles involved in the estrous cycle and FecB in the sheep oviduct, RNA sequencing was performed to analyze the transcriptomic profiles of mRNAs and long noncoding RNAs (lncRNAs) in the oviduct during the estrous cycle of STH sheep with mutant (FecBBB) and wild-type (FecB++) genotypes. In total, 21,863 lncRNAs and 43,674 mRNAs were screened, the results showed that mRNAs had significantly higher expression levels than the lncRNAs, and the expression levels of these screened transcripts were lower in the follicular phase than they were in the luteal phase. Among them, the oviductal glycoprotein gene (OVGP1) had the highest expression level. In the comparison between the follicular and luteal phases, 57 differentially expressed (DE) lncRNAs and 637 DE mRNAs were detected, including FSTL5 mRNA and LNC_016628 lncRNA. In the comparison between the FecBBB and FecB++ genotypes, 26 DE lncRNAs and 421 DE mRNAs were detected, including EEF1D mRNA and LNC_006270 lncRNA. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes functional enrichment analysis indicated that the DE mRNAs were enriched mainly in terms related to reproduction such as the tight junction, SAGA complex, ATP-binding cassette, nestin, and Hippo signaling pathway. The interaction network between DE lncRNAs and DE mRNAs indicated that LNC_018420 may be the key regulator in sheep oviduct. Together, our results can provide novel insights into the oviductal transcriptomic function against a FecB mutation background in sheep reproduction.

Ovine KRTAP36-2: A New Keratin-Associated Protein Gene Related to Variation in Wool Yield.

Keratin-associated proteins (KAPs) are structural components of wool fibres. High-glycine/tyrosine (HGT)-KAPs are a subset of the KAP family, and their abundance in fibres varies. In this study, we report the discovery of an ovine HGT-KAP gene to which we assigned the name KRTAP36-2. Polymerase chain reaction and single-strand conformation polymorphism (PCR-SSCP) analyses revealed four variants of this gene in a screening population of 170 sheep from a variety of breeds. The DNA sequencing of the variants revealed four single-nucleotide polymorphisms (SNPs) and a dinucleotide deletion. Three of these SNPs were in the coding region, and one of these was non-synonymous and potentially led to the amino acid substitution p.Cys27Gly near the middle of the protein. The remaining SNP was located near the putative TATA box, and the di-nucleotide deletion was near the putative transcription initiation site. The effect of this variation in KRTAP36-2 was investigated in 274 Southdown × Merino lambs that were the progeny of five sires. Variation was only found to be associated with wool yield, that is, the proportion of the greasy fleece that remained as clean fleece upon scouring (expressed as a percentage). This may have some value in increasing wool production.

Integrated Multi-Tissue Transcriptome Profiling Characterizes the Genetic Basis and Biomarkers Affecting Reproduction in Sheep (Ovis aries).

The heritability of litter size in sheep is low and controlled by multiple genes, but the research on its related genes is not sufficient. Here, to explore the expression pattern of multi-tissue genes in Chinese native sheep, we selected 10 tissues of the three adult ewes with the highest estimated breeding value in the early study of the prolific Xinggao sheep population. The global gene expression analysis showed that the ovary, uterus, and hypothalamus expressed the most genes. Using the Uniform Manifold Approximation and Projection (UMAP) cluster analysis, these samples were clustered into eight clusters. The functional enrichment analysis showed that the genes expressed in the spleen, uterus, and ovary were significantly enriched in the Ataxia Telangiectasia Mutated Protein (ATM) signaling pathway, and most genes in the liver, spleen, and ovary were enriched in the immune response pathway. Moreover, we focus on the expression genes of the hypothalamic-pituitary-ovarian axis (HPO) and found that 11,016 genes were co-expressed in the three tissues, and different tissues have different functions, but the oxytocin signaling pathway was widely enriched. To further explore the differences in the expression genes (DEGs) of HPO in different sheep breeds, we downloaded the transcriptome data in the public data, and the analysis of DEGs (Xinggao sheep vs. Sunite sheep in Hypothalamus, Xinggao sheep vs. Sunite sheep in Pituitary, and Xinggao sheep vs. Suffolk sheep in Ovary) revealed the neuroactive ligand-receptor interactions. In addition, the gene subsets of the transcription factors (TFs) of DEGs were identified. The results suggest that 51 TF genes and the homeobox TF may play an important role in transcriptional variation across the HPO. Altogether, our study provided the first fundamental resource to investigate the physiological functions and regulation mechanisms in sheep. This important data contributes to improving our understanding of the reproductive biology of sheep and isolating effecting molecular markers that can be used for genetic selection in sheep.